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Abstract 

Rigorous QCD predictions for decay rates of the P-wave states of lieavy quarkonia 
are presented. They are based on a new factorization theorem which is vahd to leading 
order in the heavy quark velocity and to all orders in the running coupling constant 
of QCD. The decay rates for all four P states into light hadronic or electromagnetic 
final states are expressed in terms of two phenomenological parameters, whose coeffi- 
cients are perturbatively calculable. Logarithms of the binding energy encountered in 
previous perturbative calculations of P-wave decays are factored into a phenomenolog- 
ical parameter that is related to the probability for the heavy quark-antiquark pair to 
be in a color-octet S-wave state. Applying these predictions to charmonium, we use 
measured decay rates for the Xci and Xc2 to predict the decay rates of the XcO and he- 



One of the earliest applications of perturbative QCD was the calculation of the light 
hadronic and electromagnetic decay rates of S-wave quarkonia That calculation was 
based on the assumption that the annihilation of the heavy quark and antiquark is a short 
distance process which, because of the asymptotic freedom of QCD, can be computed in 
perturbation theory. It was assumed that nonperturbative effects could be factored into 
Rns{^), the nonrelativistic wavefunction at the origin. This assumption has been supported 
by subsequent calculations beyond leading order 0, |^. Taking a similar approach in the 
case of P-wave states |^ , one might expect to be able to calculate their decay rates into light 
hadrons in terms of a single nonperturbative input R'^p{0), the derivative of the nonrela- 
tivistic wavefunction at the origin. Unfortunately explicit calculations at order (leading 
order for ^Pi and ^Pi (Ref. ^j) and next-to-leading order for ^Pq and '^P2 (Ref. 0)) reveal 
infrared divergences — a clear indication of sensitivity to nonperturbative effects beyond 
those contained in R'^p{0). In previous phenomenological applications, the divergence has 
been replaced by a logarithm of the binding energy or the confinement radius ^ or the 
radius of the bound state None of these prescriptions has any fundamental justification. 

A rigorous treatment of the P-wave decays requires a clean separation between short 
distance effects, which can be calculated as a perturbation series in the running coupling 
constant of QCD, and long distance effects. The long distance effects must either be calcu- 
lated with some nonperturbative method, such as lattice QCD, or else absorbed into a small 
set of parameters that can be determined phenomenologically. In this paper, we present 
predictions for P-wave decays based on a rigorous QCD analysis. The QCD predictions are 
based on factorization theorems |^ that will only be justified at an intuitive level in this 
paper, since our focus will be on their phenomenological implications. The decay rates of 
the four P states into light hadrons, as well as the decay rates of the ^Pq and ^P2 states into 
two photons, are all predicted in terms of two phenomenological parameters. In addition to 
R'^p{0), there is a second parameter associated with the probability for the heavy quark and 
antiquark to be in a color-octet S state. When these predictions are combined with equally 
rigorous relations between the rates for radiative transitions between P and S states, their 
predictive power is quite remarkable. From the decay rates of the ^P2 and ^Pi states of the 
first radial excitation of quarkonium, we are able to predict the inclusive decay rates of the 



^Po and ^Pi states. For charmonium, some of our predictions differ significantly from the 
presently accepted results. 

The n'th radial excitation of heavy quarkonium is split into angular-momentum states 
fi(^^+^Lj) with parity (—1)^+^ and charge conjugation (—1)^+"^. Our basic results are fac- 
torization formulae for the decay rates of these states into light hadronic or electromagnetic 
final states. The factorization formulae are valid to leading order in f ^, where v is the typical 
velocity of the heavy quark, and to all orders in the QCD coupling constant ^^(Mq). For 
S-wave (L = 0) and P-wave (L = 1) states, they have the schematic form 

r (nC'+'S) ^X) = Gi(n) f i (QQC'^'S) ^ x) , (1) 

r (n(2^+ip) ^ X) = i7iHfi(gQ(25+ip)^x) + Hs{n)fs{QQC''''S)^x) . (2) 

In ([I|) and (^, all nonperturbative effects are factored into the parameters Gi{n), Hi{n), 
and HqIu). They depend on the radial quantum number n, but are independent of the total 
spin 5* and total angular momentum J to leading order in v"^. The spin- independence follows 
from the fact that the QCD interactions of heavy quarks are independent of the spin of the 

quark, up to relativistic corrections. The factors Fi and Fg are hard subprocess rates for 
the annihilation of the heavy quark and antiquark at threshold, with the QQ in a '^^'^^L 
angular-momentum state and in a color-singlet state for Fi and a color-octet state for Fg. 
The subprocess rates are calculable as a perturbation series in q;s(Mq), the QCD running 
coupling constant evaluated at the heavy-quark mass. The parameters Gi, Hi, and Hg are 
proportional to the probabilities for the bound state to contain a QQ pair in a color-singlet 
S-wave, a color-singlet P-wave, and a color-octet S-wave state, respectively, with separation 
r — 0. They may simply be taken as phenomenological parameters to be determined by 
experiment, but they can also be given rigorous nonperturbative definitions and can 
therefore be measured using lattice simulations of QCD. In the factorization formula (^, 

and Fi depend on an arbitrary factorization scale fi in such a way that the complete decay 
rate is independent of /i. We have implicitly set fi = Mq. 

The color-octet contributions to P-wave decays {HsTg) may seem peculiar if one is 
used to thinking of mesons as color-singlet QQ states. However, any meson is a superposition 



of many components, involving any number of quarks and gluons: 

W) = i^QQ\QQ) + i'QQ.mg) + •••• (3) 

In some of these components, notably \QQg), the heavy quark and antiquark are in a color- 
octet state. The color-singlet and color-octet pieces of the our factored decay rates represent 
contributions coming from QQ and QQg components of the meson respectively. The proba- 
bility carried by the QQg state (and higher states) is of order v"^ — heavy quarks don't easily 
radiate a gluon. Consequently such components can be neglected for many applications, in- 
cluding S-wave annihilations. However, the \QQ) contribution to annihilations of P-wave 
quarkonium is suppressed by v^, owing to the angular- momentum barrier, which pushes the 
quarks apart. Furthermore, in the QQg component of P-wave quarkonium, the quark and 
antiquark can be in an S-wave state, with no angular-momentum barrier to hinder their an- 
nihilation. The QQg contribution to the annihilation therefore competes with, and in some 
cases even dominates, that coming from QQ. The enhanced role of the QQg component 
makes the decays of P states particularly interesting: they are among the very few processes 
in quarkonium physics that give us a glimpse of physics beyond the simple quark potential 
model. 

To leading order in v'^, Hi is directly related to the nonrelativistic wavefunction 

V'qq = Rnp{r)Yim{ry. 



There is no simple formula for Hg in terms of Rnp since Hg is determined by the QQg 
wavefunction ipQQg rather than by iPqq- However, in the limit of very large quark mass, 

Hg is dominated by QQg configurations in which the gluon has large momentum, and an 
approximate perturbative relation ^ can be obtained between Hi and Hg. 

where En is the binding energy, /?o = (33 — 2nf)/6, and nj is the number of light flavors. 
The constant accompanying the logarithm in (Rl) cannot be computed perturbatively and is 



probably quite important for charmonium and bottomonium. 

If in calculating the gluon emission process that gives (|^), one neglects the running 
of the coupling constant, then the perturbative expression for reduces to 

H^{n) ~ }^aM(^\ H,{n). (6) 



277r \En 

This logarithm of the binding energy is precisely the infrared divergence that was found in 
previous (nonrigorous) perturbative analyses of P-wave decays j^. Infrared divergences 

arose in this earlier work because the QQg component of the meson was neglected. In our 
analysis the infrared sensitivity is factored into the nonperturbative parameter ifg, so the 

subprocess rate Fi involves only hard contributions. 

For the decay rates of the P states into light hadrons, the factorization formula (0) 
in more explicit form is 

T(n{^Vj)^l.h) = H^{n)f^{QQi^FJ)^partons) 

+ Hs{n) fg (gg('Si) ^ partons) , J = 0, 1,2 , (7) 



F (n(iPi) ^ l.h.) = Hiin) fi (gg(^Pi) ^ partons 

+ H^{n) fg [QQC^o) partons 



where "/./i." on the left side of (|^) or (|^) represents all final states consisting of light hadrons 
and ^^partons" on the right side represents perturbative final states consisting of gluons and 
light quark-antiquark pairs. Note that, for the ^Pj decays, the second term on the right side 
of d^) is independent of J. For the decay rates of the ^Pq and ^P2 states into two photons, the 
factorization formula (0) simplifies to the color-singlet term only, because color conservation 
forbids the annihilation of the QQ in a color-octet state: 

F (nCPj) ^ 77) = H,{n) f, (QQCPj) ^ 77) , ^ = 0, 2 . (9) 



When applied to decays into a hard photon plus light hadrons, the factorization 
formula (§) has remarkable implications. The color-octet term allows the ^Pi state to decay 



into a hard photon plus hght hadrons at order ^^^(Mq) through the subprocess QQ(}So) 
'yg. This decay produces a final state consisting of a hard photon recoiling against a hard 
gluon jet and against the hadrons from the fragmentation of the gluon in the QQg state. It 
is possible that this dramatic decay mode of the ^Pi state could serve as an experimental 
signature for this particle. 

The leading subprocess rates Fi and Fg for light-hadronic and electromagnetic decays 
of P states are listed in Table 1. They can be extracted from previous calculations of P-wave 
decays 0, by using the expressions (0) and for Hi and Hg. Since next-to- leading- 
order corrections in have not been computed for the Fg's, we work only to leading order 
throughout. However, we note that, since Hi and Hg are independent parameters, a higher- 
order contribution involving Hi could, in principle, be numerically important compared to 
a leading-order contribution involving Hg. 

The two nonperturbative parameters Hi and Hg can be obtained by measuring the 
decay rates into light hadrons of the ^Pi and ^P2 states, which are the two P states most 
accessible to experiment. At leading order in and q;s(Mq), they are 

^ 45 F (3P, ^ /./..) -F(3P,^ /./..) 

Hs . -^'-^^ . (11) 

For applications at leading order in ^^(Mq), it is convenient to eliminate Hi and 
Hs to obtain direct relations between decay rates. From (|^ and (^, we get the following 
relations between the decay rates of the spin-triplet states, valid to leading order in f ^ and 
in ^^(Mq): 



These differ from previous predictions 0], wliicli can be obtained by setting r('^Pi l.h.) = 
0. At leading order, the decays of the two spin-1 states involve only the color-octet term, so 
the ratios of the decay rates are simply the ratios of the subprocess rates Fg in Table 1: 



r(3Pi ^ l.h.) In 



f 



The prediction (|15D was first made by Barbieri et al. 

There are also relations among the radiative transitions between P and S states that 
are accurate up to relativistic corrections of order v"^ (Ref. 0). The predictions for the decay 
rates of P states into S states plus a photon are 

reP^^r(3p,^_^^^_^_^_ (IT) 

7 7 

where in the denominator is understood to be the energy of the photon for the transition 
in the numerator. In terms of the masses Mp and Ms of the bound states, = {Mp — 
Ml)/{2Mp). 

We now apply our QCD predictions to the charmonium system. For the first radial 
excitation, the ^Pj states are called Xcj and the ^Pi state is called the he. We use measured 
decay rates of the Xci and Xc2 to predict the inclusive decay rates of the Xco and he. It is 
important to have reasonable estimates for the theoretical errors in our predictions if they 
are to be compared with experimental results. The two main sources of theoretical error 
are relativistic corrections and higher-order perturbative corrections. In potential models 
of charmonium, the average value of v"^ is found to be about 0.23 (Ref. JlOl)- Since our 
factorization formulae are valid only to leading order in v"^, we expect an error on the order 
of 20% due to relativistic effects. Similarly, we expect deviations on the order of 20% from 
the equalities (|l^) involving the radiative transition rates. We can estimate the perturbative 
error from the size of the perturbative corrections in other bound state calculations. Based 
on a number of next-to-leading-order calculations for S-wave bound states we estimate 



the perturbative error to be 4Q;s(Mc)/7r, where Mc is the mass of the charm quark. To avoid 
the ambiguity in the value of Mc, we determine ^^(Mc) by taking the couphng constant 
Q;s(Mfe) = 0.179 ± 0.009 extracted from bottomonium decays [0 and evolving it down to the 
scale Mf.. This does not require that we know and Mf, separately, but only their ratio, 
for which we use the ratio of the J/ip and T masses: Mc/M^ ~ 0.33. The resulting value 
of the coupling constant is as{Mc) = 0.25 ± 0.02. Our estimate of the perturbative error is 
therefore 30%. We treat the 8% error in the value of as{Mc) itself as an experimental error. 
We take the QED coupling at the scale Mc to be a = 1/133.3. 

We assume in our analysis that the decays into light hadrons and the radiative tran- 
sitions to J/ip or rjc are the only decay modes which contribute appreciably to the total 
decay rates of Xcj and he- In particular, we neglect pionic transitions of the P states to the 
S states, of which the most important decay modes should be J/ip + tttt and rjc + vrvr. The 
rate for the particular decay he —>■ J/ip + tttt has been estimated within a well-developed 
phenomenological framework |TI[] to be on the order of 6 keV, and the other two-pion tran- 
sition rates should be of the same order of magnitude. The errors due to neglecting these 
contributions to the total decay rates are negligible compared to other errors. 

Precision measurements of the total decay rates of the '^Pi state Xci and the ^P2 state 
Xc2 have recently been carried out at Fermilab by the E760 collaboration |]T2|. Their results, 
with statistical and systematic errors added in quadrature, are listed as input values in the 
first column of Table 2. Previous experiments have measured the branching fractions for the 
radiative transitions of the Xci and Xc2 into the J/ ip , and they are listed as input values 
in the first column of Table 3. We use the radiative branching fractions and the total decay 
rates to obtain the partial rates given in Table 2 for light-hadronic and radiative decays of 
the Xci and Xc2- 



Inserting the partial rates into light hadrons into (|10D and ([TT|), we determine the 
nonperturbative parameters for P-wave decays of charmonium to be 

Hi ~ 15.3 ± 3.7 MeV (±36%) , (18) 

Hs ~ 3.26 ± 0.73 MeV (±36%) . (19) 
Here, and throughout the remainder of this paper, the first error is from the uncertainties 



in our experimental inputs and the second error is our estimate of the theoretical uncer- 
tainty. The theoretical error is computed by combining the relativistic error of 20% and the 
perturbative error of 30%, using the standard formulae for propagating independent errors. 
The ratio of the two nonperturbative parameters in (0) and (|19D is H^/Hi 0.21. This is 
roughly consistent with the value that one would obtain from the leading-log perturbative 
expression (^) by arbitrarily setting as{En) ~ 1. However, there is no a pnon justification 
for using the leading-log approximation at such small values of the heavy quark-mass. 

We proceed to calculate the total decay rates of the Xco and he- The experimental 
values for the radiative transition rates of Xci and Xc2 in Table 2 are in good agreement with 
the theoretical prediction (|T^. This gives us confidence in our predictions in Table 2 for the 
radiative transition rates of Xco and he- In calculating the radiative transition rate of the 
he, we have assumed that its mass is at the center of gravity of the nine spin states of Xco, 
Xci, and Xc2- = 3525 MeV. The predictions in Table 2 for the decay rates of Xco and 
he into light hadrons are obtained by using (|12|) and (p^Sf). Adding the light-hadronic and 
radiative decay rates, we obtain the total decay rates given in Table 2. The prediction for 
the Xco differs from the presently accepted value [|13] of (14 ± 5) MeV by several standard 
deviations. Our calculations suggest a real discrepancy between this measurement of the Xco 
decay rate and the experimental data on the Xci and Xc2 that we have used as our input. 

In Table 3, we present our predictions for the branching fractions for decays of the 
charmonium P states into final states containing photons. The branching fractions for Xco 
and Xc2 into two photons are calculated from (|TB|) and (0). The branching fraction for the 
he to decay into a hard photon plus light hadrons is calculated using (|T^) . Our prediction 
for the radiative branching fraction of the Xco is significantly larger than the accepted value 



[|T^ of 0.0066 ±0.0018. Our prediction for the branching fraction of the Xc2 into two photons 
is considerably smaller than the accepted value |jl3| of (11 ± 6) x 10^"^. A rather tantalizing 
prediction is that the branching fraction for the he to decay into a hard photon plus light 
hadrons is about 2%. This is large enough that it may be possible to detect the decay of the 
he by observing the resulting hard photon. 

P-wave annihilations provide unique information on the dynamical role of the gluon in 
determining hadron structure. They challenge us to go beyond the quark-potential picture in 



modelling hadrons. In this paper we have outlined the first rigorous formalism for describing 
these processes in QCD. We have applied the formalism in a detailed analysis of charmonium 
P-wave decays. Within our formalism it is possible to refine systematically the theoretical 
predictions for the decay rates, both by computing higher-order perturbative contributions, 
and also by using lattice simulations to compute the nonperturbative parameters. The 
promise of significant improvements in both theory and experiment make the P-wave decays 
important testing grounds for ideas about both perturbative and nonperturbative QCD. 
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Table Captions 

1. Subprocess rates Fi and Fg for decays of P-wave quarkonium states. 

2. Predictions for total and partial decay rates of P-wave charmonium states. 

3. Predictions for branching fractions of P-wave charmonium states. 
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r in MeV 


r{i.h.) 


r(7j/V'),r(7^c) 


XcO 


4.8 ± 0.7 (±35%) 


4.7 ±0.7 (±36%) 


0.099 ±0.010 (±20%) 


Xcl 


INPUT: 0.88 ±0.14 


0.64 ±0.10 


0.240 ± 0.041 


Xc2 


INPUT: 1.98 ± 0.18 


1.71 ±0.16 


0.267 ±0.033 


he 


0.98 ± 0.09 (±22%) 


0.53 ± 0.08 (±36%) 


0.45 ± 0.05 (±20%) 



Table 2 







B{^+l.h.) 


5(77) 


XcO 


0.021 ± 0.004 (±40%) 




(6.8 ±1.9) X 10-^ (±50%) 


Xcl 


INPUT: 0.273 ±0.016 







Xc2 


INPUT: 0.135 ±0.011 




(4.1 ± 1.1) X 10-4 (±36%) 


he 


0.46 ± 0.05 (±22%) 


0.017 ±0.003 (±42%) 
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